Gate-induced interlayer asymmetry in ABA-stacked trilayer graphene 
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We calculate the electronic band structure of ABA-stacked trilayer graphene in the presence of 
external gates, using a self-consistent Hartree approximation to take account of screening. In the 
absence of a gate potential, there are separate pairs of linear and parabolic bands at low energy. 
A gate field perpendicular to the layers breaks mirror reflection symmetry with respect to the 
central layer and hybridizes the linear and parabolic low-energy bands, leaving a chiral Hamiltonian 
essentially different from that of monolayer or bilayer graphene. Using the self-consistent Born 
approximation, we find that the density of states and the minimal conductivity in the presence of 
disorder generally increase as the gate field increases, in sharp contrast with bilayer graphene. 

PACS numbers: 71.20.-b,81.05.Uw,73.63.-b,73.43.Cd. 



Pioneering experiments [ll, [3, H, |3| demonstrated 
graphene-based transistors using a back gate to vary the 
carrier density continuously from electron to hole chan- 
nels, with a minimal conductivity for nominally-zero car- 
rier density. The switching of a graphene-based transistor 
would be improved by opening an energy gap between 
the conduction and valence bands, possibly bv lateral 
confinement of electrons in etched structures [a, @i Hi [1| 
or by employing gates to induce interlayer asy mmetry in 
bilayer graphene [1, [IS [ll|, [11 [ll, [ll lli, [3 . Recently, 
experimental attention has turned towards the properties 
of ABA-stacked trilayer graphene, Fig.^a), [l3,[ll[ll. 
Theory suggests that the bands are of two separate types 
[H [nj, [la, [m, [H, [Hi: two almost-linear bands rem- 
iniscent of the bands in monolayer graphene and four 
parabolic bands similar to those in bilayer graphene. 
This raises the expectation that the electronic behavior 
will display no new features as compared to monolayer 
or bilayer graphene. 

In this paper, we show theoretically that the response 
of ABA-stacked trilayer graphene to external gate poten- 
tials is in fact qualitatively different from that in mono- 
or bi-layer graphene. We use an effective-mass model to 
self-consistently determine the electronic band structure 
and we show how the breaking of mirror reflection sym- 
metry by interlayer asymmetry causes hybridization of 
the linear and parabolic bands. Rather than opening a 
gap, as in bilayer graphene this leaves two bands near 
zero energy which support chiral quasiparticles. Employ- 
ing a self-consistent Born approximation to estimate the 
minimal conductivity as a function of interlayer asym- 
metry, we find that the conductivity generally increases 
as asymmetry increases, in sharp contrast with bilayer 
graphene as illustrated in Fig. (Hb) . 

A description of trilayer graphene in the presence of ex- 
ternal gates must include two parameters that take into 
account differences in the potentials Vi, V2, and V3 of the 
three layers. The first, Ai = —e{Vi — V3)/2, describes 
the average energy difference between each adjacent layer 
[l3,[il|,lll, while the second, A2 = -e{Vi - 2V2 + Va)/^, 
describes the difference between the energy of the cen- 
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FIG. 1: (a) Schematic of the ABA-stacked trilayer lattice 
containing six sites in the unit cell, A (white circles) and B 
(black circles) on each layer, showing the Slonczewski-Weiss- 
McClure parameterization [24] of relevant couplings 70 to 75. 

(b) The conductivity versus external asymmetry A^'^"*'' , calcu- 
lated for trilayer and bilayer graphene using the self-consistent 
Born approximation and the band model including 70 and 71 . 

(c) Schematic of trilayer graphene (three thin black lines at 
X — —d, 0, d) with top and bottom gates (thick black lines at 
X = Lt, —Li,) separated from the trilayer by dielectric media 
(gray shaded areas). 



tral layer and the average of the outer layers. We model 
the effect of back and top gates by considering the tri- 
layer as three conducting parallel plates as illustrated in 
Fig- [He), with respective electron densities ni, 712 and n^, 
located at x = —d, 0, and -\-d, respectively, where d is 
the interlayer spacing, and the permittivity of the trilayer 
interlayer spaces (without the screening effect of 7r-band 
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electrons of the trilayer graphene) is Sr- The back (top) 
gate at X = —if, {x = +Lt), held at potential Vb (Vt), is 
separated from the trilayer by a dielectric medium with 
relative permittivity £{, (st)- Using elementary electro- 
statics, we relate the external gate potentials, the elec- 
tron densities on the layers, and the interlayer asymmetry 
parameters: 



£bVfc 
Lb 

Ai = 



A, 



StVt , 
—— = e [ni 



"2 + ns) 



StVt 
Lt 

6er 



£bVb 



+ ^("1 -"3), 
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In the following, we use the total electron density ntot = 
ni -1-^2-1-713 and A^'^''*^ as external parameters instead of 
Vt and Vg, where A^'^'^'^ = StVt/Lt — SbVb/Lb is the value 
of Ai that would occur if screening were negligible. 

We model ABA-stacked trilayer graphene as three cou- 
pled honeycomb lattices including pairs of inequivalent 
sites {^1,51}, {A2,B2}, and {A3,B3} in the bottom, 
center, and top layers, respectively. The layers are ar- 
ranged according to Bernal (A-B) stacking |23|, Fig.[T][a), 
such that sites Bl, A2, and B3 lie directly above or 
below each other. We employ an effective-mass model 
adopting the Slonczewski-Weiss-McClure parameteriza- 
tion, [2J] where each parameter is related to relevant 
coupling in the tight-binding model: 70 describes nearest- 
neighbor {Ai-Bi for i = {1,2,3}) coupling within each 
layer, 71 describes strong nearest-layer coupling between 
sites {B1-A2 and A2-B3) that lie directly above or below 
each other, 73 (74) describes weaker nearest-layer cou- 
pling between sites A1-B2 and B2-A3 {A1-A2, B1-B2, 
A2-A3, and B2-B3). With only these couplings, there 
would be a degeneracy point at each of two inequivalent 
corners, K±, of the hexagonal Brillouin zone ^25^ but 
this degeneracy is broken by next-nearest-layer coupling 
72 (between AI and A3), 75 (between Bl and -63) and S, 
which is the on-site energy difference between AI, B2, A3 
and Bl, A2, B3. Note that the parameter A often used 
in models of three-dimensional (3D) graphite is given by 
A = (5 -|- 72 — 75. In trilayer graphene, the presence of 
a surface may induce a modification in the value of the 
band parameters as compared to those in bulk graphite. 
Here, parameter A2 = —e{Vi — 2V2 + V3)/Q, takes into 
account a possible difference between the energy of the 
central layer and the average of the outer layers, and, 
in general, surface effects may contribute to a non-zero 
value of A2 . 

In a basis with components ipAi, '4'Bi, ipA2 , "0-82, ''Pas, 
ipB3, the ABA-stacked trilayer Hamiltonian is 
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where operator tt = ^p^ + ip ^, is related to the in- 
plane momentum p — {px,Py) [25|, effective velocities 
are v = [■\/3/2)a'yo/h, U3 — {^/3/2)a'yz/K and W4 = 
(^3/2)074/?!, Ui = ~eVi, and ^ = ±1 is the valley in- 
dex K±. Exploiting mirror reflection symmetry of the 
lattice in the plane of its central layer, Fig.[TJa), we per- 
form a unitary transformation to a basis consisting of 
linear combinations of the atomic orbitals [22i |. namely 

[V'Al - ^A3]/\/2, [V'Sl - ■(/'-B3]/\/2, [V'Al + V'A3]/a/2, 1pB2, 
1pA2, bpBl + tl^B3]/V2: 
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where the average on-site energy [Ux + U2 + C^3]/3 has 
been set equal to zero. The Hamiltonian H has a 2 x 2 
block Hjn and a 4 x 4 block Hb on the diagonal, connected 
by a simple off-diagonal block D. Block H,n is similar to 
the Dirac-type Hamiltonian of monolayer graphene and it 
contributes two bands near zero energy whereas block Hb 
is reminiscent of the Hamiltonian of bilayer graphene 0] , 
except that terms proportional to 71, 73, and 74 appear 
with a factor -\/2 22]. The latter gives two bands split 
away from zero by energy ±-\/27i and two bands near 
zero energy. 

The monolayer-like block has wave functions possess- 
ing odd mirror reflection symmetry, while the wave func- 
tions of the bilayer part are even. Since the interlayer 
asymmetry Ai is the only parameter that breaks mir- 
ror reflection symmetry, its role is qualitatively different 
from the other parameters, coupling the monolayer-like 
and bilayer-like blocks. For large Ai, two of the low- 
energy bands, related to orbitals [ipAi — ^Az\l^ and 
[V'Ai + ■0^3] /v^, split away from zero by energy e « ±Ai 
at the K point, leaving only two bands near zero, associ- 

ated with ^f' = — '(/'-B3]/'s/2, '0B2) ■ To obtain an 

approximate Hamiltonian H^s for 4"', we denote H2 as 
the diagonal block of Hamiltonian H corresponding to 
these two low-energy components, H,^ as the 4x4 diago- 
nal block corresponding to the high-energy components, 
and V as the off-diagonal 2x4 block coupling H2 and 
i?4. The Schrodinger equation for can be expanded 
up to first order in e as \H2 - VHl^V'^\^' = eS**' with 
S" = 1 -I- VH^'^V'^ . Then, the effective Hamiltonian for 
* = S^l"^^' becomes H^s ~ S-^l'^\H2-VHl^V'^\S-^l'^ . 
For the moment, we focus on the role of Ai by consider- 
ing A2 = 72 = 73 = 74 = 75 = = 0. For large enough 
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Ai (|7i| > |Ai| > |e|), iJcff is written 



as 



cff 



X 





X 

V271 
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For plane-wave eigenstates at zero magnetic field, tttt^ 
is just a number, p^. The first factor of tt = £,px + ipy 
in operator X ensures that such eigenstates are chiral. 



/\/2 with = tan ^{py/px) 



The expression for the eigenenergies is 
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which generalizes Eq. (22) of Ref. [ll|, showing that there 
is a small overlap (5e ~ Af/71 between the two low- 
energy bands that cross at p ~ Ai/v ^1^. This behavior 
contrasts with that of bilayer graphene, where interlayer 
asymmetry introduces an energy gap between the low- 
energy bands Q and tends to suppress the chiral nature 
of quasiparticles in them. 

For given parameters Ai, A2 and fixed total density 
ntot, the electron densities ni,n2 and may be deter- 
mined by summing ( | TpAi P + 1 V'Bi P ) /^^ over the occupied 
eigenstates of the Hamiltonian However, such densi- 
ties are also related to Ai and A2 through Eqs. ^ and 
([3]), so it is necessary to solve this set of equations self- 
consistently in order to obtain values of Ai and A2 for 
given external parameters ntot and A^''''*\ A similar pro- 
cedure has been applied to bilayer graphene [H, [H, 
and to many-layered graphene 26] . This Hartree approx- 
imation neglects effects including exchange interaction, 
possible deformation of atomic orbitals in the applied 
electric field, and the role of a orbitals in screening, but 
comparison with density functional theory ^13j in bilayers 
suggests that it is qualitatively accurate. 

For 72 = 73 = 74 = 75 = i5 = with ntot = 0, it 
is possible to perform a linear response calculation for 
infinitely small A''''^*^ Within the first order in Ai, we 
have ni — 77.3 = IIAi with 
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where gs = gv — 2 are the spin and valley degeneracies, 
respectively, \a) and the eigenstates and eigenen- 

ergy of the Hamiltonian without Ai or A2, /(e) is the 
Fermi distribution function with zero Fermi energy. Us- 
ing Eqs. ^ and ([3]), we obtain the self-consistent solution 
Al ~ A^'^'^^^ecff with eoff = 1 — (e^(i/2ej.)H. A2 is never 
induced. Typical parameters v = 1.0 x 10^ m/s, 71 = 0.4 
eV, d = 0.334 nm, e^. = 2 give l/eog ~ 0.61. 

To determine the band structure taking into account 
all the parameters, we find Ai and A2 self-consistently 
by employing an iterative numerical approach. We first 
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FIG. 2: Self-consistently calculated band structures in tri- 
layer graphene near the K point, with ntot = 0. Left plots 
are for the model including only 70 and 71, with (a) no 
asymmetry A^'^'^*' — and (b) finite external asymmetry 
^(cxt) _ g Right plots are for the full parameter model 
including 7^(1 = 1,2, ■••5), 5, with (c) A^'^"*' = and (d) 
0.571. Dashed and solid curves represent 6 = Q and 7r/6. The 
self-consistently calculated value of Ai is shown in the lower 
side of each plot. The thin horizontal line shows the Fermi 
energy 



use Al = aJ^'^*'' and A2 = as initial values in the 
Hamiltonian Eq. ([4]) and determine the Fermi energy so 
that the total density is equal to ntot ■ Then we calculate 
Ui (i = 1, 2, 3) from the occupied eigenstates, which give 
a new set of Ai and A2 through Eqs. ^ and We 
iterate this process until Ai and A2 converge. 

Figure [5] (a) and (c) show the self-consistent band 
structures at zero external field A*^''''*) — and zero 
doping ntot = 0. To illustrate the role of the extra 
band parameters we compare (a) the simple model in- 
cluding only 7o, 71 and (c) the full-parameter model with 
72 = -O.O571, 75 = O.I71, (5 = 0.12571, U3(cx 73) = O.li) 
and «4((x 74) = 0.0141; (typical values quoted for bulk 
graphite [23|). The plots show the vicinity of zero energy, 
covering the monolayer-like band and the lower branches 
of the bilayer- like band. In (c), we see that 72, 75 and 
5 shift the center of the monolayer-like band upward in 
energy relatively to the bilayer-like band. Also, the trig- 
onal warping effect due to 73 is observed as a difference 
between 6 = tan~^(pj,/pa:) = and 7r/6 [9*]. Figures [2jb) 
and lUd) display the corresponding plots in the presence 



of a finite external field A.\ 



(oxt) _ 



0.571. The values of Ai 



determined self-consistently are shown in the lower side 
of each plot. In every case the screening ratio Ai/A^'''^*'' 
is about 0.6, which is close to the linear response theory. 
For Fig. [2] (b), where only 70 and 71 are included, there 
is a small overlap at zero energy described by Eq. ([H]) . In 
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FIG. 3: (a),(c) Density of states and (b),(d) conductivity at 
J^tot = as functions of A^™'' for the model including only 
7o and 71 (left) and the full parameter model (right). Solid, 
dotted, and dashed lines are for progressively larger disorder 
strengths. 



the full parameter model (d), there is a similar amount 
of band overlap while the exact magnitude of momentum 
at the crossing point vp ~ Ai varies with angle in a 
trigonal manner, and there is a tiny gap at those crossing 
points. In (d), the self-consistent calculations yield tiny 
A2 < O.OI71 due to non-zero 72, 75 and S. 

For each band structure we estimate the density of 
states (DOS) and conductivity using the self-consistent 
Born approximation [23, HI] . We assume that the scat- 
terers are on-site potentials localized on each layer, which 
is modeled by ^(r) = X]m=i 2 3 Si UiS{r—ri)P'^"^^ where 
Ui and Ti = {xi,yi) are the amplitude and the two- 
dimensional position of the i-th scatterer, respectively, 
and is the projection operator onto the m-th layer. 
We neglect inter- valley scattering between K±. The dis- 
order strength is characterized by = rii^pU^ / [A-kH'^v'^) 
[13] where nimp is the total number of scatterers over all 
the layers, and u — {uf). For the model with only 70 
and 71, the energy scale for the level broadening at zero 
energy is given by r ~ {tt/V2)Wji. Following Ref. 113, 
we compute the self-energy and the vertex corrections 
for the velocity operators, and calculate the conductivity 
using the Kubo formula. 

Figures ^a,) and ^c) show the DOS as a function 

of Af'^'^ at ntot = 0, for several values of the dis- 
order strength W. The left and right panels corre- 
spond to the simple model with 70 and 71, and the 
full parameter model discussed previously, although the 
behavior in each case is similar, DOS increases with 
^^cxt)^ Figure [3]Jb) , (d) are plots of the conductivity cor- 
responding to (a),(c) respectively. The general trend is 



for an increase of conductivity as A^'^'^*^ increases, ex- 



cept for the vicinity of A^"'^*-' = in panel (d). This 
may be roughly understood by considering the relation 
a — e^pFVpT/2 with velocity vp, DOS pp and relax- 
ation time T at the Fermi energy. When we assume that 
all the states on the Fermi energy are equally mixed by 
disorder, we have r (x Pp^, suggesting that the conduc- 
tivity is determined by Vp. The dispersion, Eq. ([5]), ap- 
proximates, in regions far from the origin vp ^ |Ai|, to 
e w ±('i;V _ 3A^/2)/(V27i), showing that the electron 
and hole bands are pushed towards zero energy by the in- 
troduction of Ai. This leads to an increase in the expec- 
tation value of the band velocity in a disorder-broadened 
energy window near e = 0, and thus the conductivity at 
the charge neutral point is enhanced. 

In the simple 70-71 model. Fig. [31Jb), the conductivity 
at A^''''*'' = takes a universal value a — ig^gsC^/ (27r^fi,) 
independently of W . This is because the Fermi energy 
coincides with the degeneracy point of the monolayer 
and bilayer bands, and the value is indeed equal to the 
summation of the minimum conductivity of monolayer 
graphene [13| and that of bilayer graphene [1^ estimated 
in the self-consistent Born approximation. In Fig. ^d) , 
for the full parameter model, the conductivity is largely 
enhanced around A^'^'^*'' — 0, because, as observed in 
Fig. m^c), the Fermi energy crosses the off-center part 
of the monolayer- like band making a large contribution 
to the typical band velocity. The conductivity drops 
sharply as A^°'^*'' grows from zero as the monolayer band 

is gapped away. When A^"'^*'' is increased further, the 
conductivity grows similarly to, but a little more slowly 
than, Fig. [3fb), because of the tiny gap at the band 
crossing point observed in Fig.[2Jd). The differences be- 
tween Fig. ^h) and (d) become smaller for larger W, 
as disorder-broadening masks details dependent on the 
precise values of band parameters. 

To conclude, we have shown that the breaking of mirror 
reflection symmetry by inter layer asymmetry Ai in ABA- 
stacked trilayer graphene causes hybridization of the lin- 
ear and parabolic bands, leaving just two bands in the 
vicinity of zero energy. The band hybridization produces 
an increase in density of states and typical band velocity 
with asymmetry Ai, leading to an increase in minimal 
conductivity in qualitative agreement with recent trans- 
port experiments [l^. As demonstrated in Fig. [TJb), 
which compares the conductivities of trilayer and bilayer 
graphene, the response of trilayers to gate-induced asym- 
metry is in sharp contrast with bilayers, where the con- 
ductivity is suppressed by a perpendicular electric field 
owing to the opening of a gap between the electron and 
hole bands . 
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